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Abstract 

The Guard-Guardee model for plant immunity describes how resistance proteins (guards) in host cells 
monitor host target proteins (guardees) that are manipulated by pathogen effector proteins. A recently sug- 



^1 gested extension of this model includes decoys, which are duplicated copies of guardee proteins, and which 

have the sole function to attract the effector and, when modified by the effector, trigger the plant immune 
response. Here we present a proof-of-principle model for the functioning of decoys in plant immunity, 
quantitatively developing this experimentally-derived concept. Our model links the basic cellular chemistry 
n I to the outcomes of pathogen infection and resulting fitness costs for the host. In particular, the model allows 



identification of conditions under which it is optimal for decoys to act as triggers for the plant immune 
response, and of conditions under which it is optimal for decoys to act as sinks that bind the pathogen 
effectors but do not trigger an immune response. 

Keywords: Plant Immunity, R protein, guard-guardee model, decoy 
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I. INTRODUCTION 

The Guard-Guardee model of pathogen recognition has become an established paradigm in 
plant immunity 121 |4J, a, [lOD. According to this model, an effector molecule, injected by the 
pathogen into a host cell to facilitate infection, targets some molecule that is important for cellular 
life. The state of this guardee molecule is monitored by a specific immune resistance (R) protein, 
also called a guard. When the state of the guardee molecule is modified by the effector, the 
R protein detects this change and activates an immune reaction that usually kills the cell, thus 
preventing pathogen proliferation and warning other cells of the threat of infection. However, 
based on recently acquired empirical data, it has been conjectured, that the plant innate immune 
system employs even more sophisticated mechanisms to fight off pathogen attacks. 

In particular, [llOl] have suggested that a duplication of some of the effector target genes and 
subsequent independent evolution of the duplicated genes can create decoys of the effector targets. 
Such a decoy is free from the non-immune cellular functions of the original target proteins and 
can evolve to become an effective bait for the pathogen effector, as well as a better trigger for its 
corresponding guard R protein. Meanwhile, in the presence of decoys the target protein might 
have fewer constraints to mutate away from the recognizability by the effector, which might be 
beneficial to the individuals without functioning R genes. Four examples of possible effector 
perception mechanisms, summarized in [10], support the decoy model. The diversity of these 
examples, which include the Pto protein as a candidate for the decoy in tomato, pBS3 in pepper, 
RCR3 in tomato as a bait for fungus effector, and the RPS2 protein in Arabidopsis, indicate that 
such decoys may often evolve independently and may be active in a wide variety of plant-pathogen 
interactions. 

It is interesting that the notion of decoy is also being used in the different, but related context of 
"decoy receptors", which are inactive duplicates of normal receptors for general signalling path- 
ways. Such decoys have high and selective binding affinity to specific ligands but lack the signal 
transduction domains of normal receptors. Decoy receptors act as sinks for ligands and prevent 
them from binding to their active targets (see, for example, [lll,l3|] or work as parts of multiprotein 



receptor complexes. 
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Importantly, it has been reported that certain cancer cells evolve to escape 



immune-cytotoxic attack by over-expressing a decoy receptor that sequesters the respective apop- 



80 



tosis ligand IJSI]. For immune reactions in plants, it therefore seems natural to assume that decoys 
can function in two different ways: first, as being baits for triggering the guard-mediated immune 



response and, second, in case of sufficient abundance and not too strong recognition by guards, by 
simply absorbing the pathogen effectors and preventing them from attacking their normal targets, 
thus preventing both infection and the immune response. In the following we will often refer to 
these two mechanisms and decoy roles as a "trigger" and "sink". 

Our goal here is to quantify the presently only qualitative description of both trigger and sink 
mechanisms of decoy functioning in host immune response, and to understand the conditions that 
might favour one or the other mechanism evolutionarily. The advantages and disadvantages of 
each of the mechanisms are summarized in the form of the overall fitness decrease of the organism 
caused by the pathogen infection. We develop a simple model which links some basic cellular 
chemistry to probabilities of possible outcomes of parasite attack and their fitness costs for a host 
organism. In our model we consider reversible binding interaction between four major types of 
molecules: pathogen effector, its target in a plant cell (guardee), the decoy of this target, and the 
guard R protein. The effector can bind both to the target and the decoy, and the guard protein 
can interact with the effector-target and effector-decoy complexes and thus trigger the immune 
response. Fitness costs are associated with the maintenance of decoy proteins, with cell death 
caused by the activated immune response, and with undetected and therefore successful infection, 
which also causes cell death, but enables pathogen multiplication and subsequent attack on a 
number of adjacent cells. We find that decoys are typically favoured evolutionarily, but depending 
on the parameters of the model, the decoy can indeed exhibit either the role of a sink for effector 
molecules, or the role of an effective trigger of the immune response. In particular, when the cost 
of maintenance of the decoy is low, higher fitness is achieved when decoys act as a sink: the decoy 
protein tends to be more abundant and simply sequester the effector without being recognized by 
the guard, thus avoiding cell death. However, when the cost of the maintenance of the decoy is 
elevated, higher fitness is achieved with the fewer decoy molecules being triggers of the immune 
response, and with a high recognizability of the effector-decoy complex by the guard protein. 

The structure of the paper is the following: In the next section we describe our model, start- 
ing with the basic cellular chemistry, linking the chemical concentrations to the probabilities of 
outcomes for the infected plant, evaluating the fitness costs of each outcome, and determining 
parameters and scenarios that maximize the fitness. The discussion and conclusion sections com- 
pletes the paper. 



II. THE MODEL 



A. Binding equilibrium 



We consider biochemical reactions in a system initially consisting of 4 types of monomers 
(proteins): pathogen effector, denoted by E, effector target (guardee) in a host cell, denoted by T, 
decoy of the effector target denoted by D, and a guard R protein denoted by R, which monitors 
the state of the effector target and its decoy. These monomers can reversibly react in the following 
ways: The effector can bind both to its target to form a complex ET and to the decoy to form a 
complex ED; the guard protein can recognize modified effector targets and modified decoys by 
binding to ET and ED complexes, thus forming RET and RED three-molecular aggregates. It 
is assumed that guard proteins do not bind to non-modified targets or decoys, and hence that the 
formation of complexes RE and RD is impossible. We assume that the chemical equilibrium 
between these binding states establishes itself fast and that only the steady state concentration of 
complexes are biologically relevant. We also assume that the concentrations of all reactants are 
uniform and therefore disregard any spatial dependence. The equilibrium concentrations of the 
constituent proteins and their complexes are uniquely defined by the law of Mass Action, ([T]), 

[E][T] = Ket[ET] (1) 

[E][D] = Ked[ED] 
[R][ET] = KretIRET] 
[R][ED] = KnEDiRED]. 

Here [A] denotes the concentration of substance A, and Kab is a dissociation constant of a com- 
plex AB. A dissociation constant has the dimensionality of a concentration, and a smaller disso- 
ciation constant means stronger binding. In addition, there are conservation laws which state that 
the total amount of each of the four proteins, i.e. the amount of free and bound protein, remains 
constant and is equal to [E]q, [T]q, [D]q, and [R]o, respectively. We assume that the values of the 
total concentrations change very slowly on the timescale of the development and detection of a 



pathogen attack and therefore approximate them by constant values. 



[E]o = [E] + [FT] + [ED] + [RET] + [RED] (2) 

[T]o = [T] + [ET] + [RET] 
[D]o = [D] + [ED] + [RED] 
[R]o = [R] + [RET] + [RED]. 



The system of 8 equations resulting from (1 1121) in principle allows to determine all 8 equilibrium 
concentrations, but the equations are nonlinear and can be solved only numerically. A simple way 
to numerically find the unique equilibrium concentrations [9] is to integrate the actual kinetic equa- 
tions that lead to the establishment of the equilibrium over time, keeping the ratio of the kinetic 
coefficients for association and dissociation equal to the corresponding dissociation constants ||7|]. 
For example, for a complex AB formed as a result of the reaction 

A + B^AB, (3) 

with a dissociation constant Kab, the corresponding system of kinetic equations is 

d[A] 



dt 

d[B] 

dt " 

d[AB] 

dt 



-[A][B]+Kab[AB] 

-[A][B]+Kab[AB] (4) 

--[A][B]-Kab[AB]. 



For each protein, the corresponding equations are numerically evolved in time until the steady 
state concentrations are reached. 

B, Probability of effector and guard success and concentration of complexes 

The steady state concentrations of protein complexes are directly related to probabilities of 
success or failure of the pathogen effector to induce changes in its target or decoy, and to the 
probability that a guard protein recognizes these changes. We assume that the extent to which the 
effector modifies the cellular environment to the benefit of the parasite depends on the fraction 
of target molecules that it binds to. Similarly, the efficiency of recognition of the effector attack 
by the guard protein depends on the fraction of guard protein molecules bound to the effector- 
target complex. While generally the dependencies between the probabilities of such effects and 



concentrations ratios can have a complex functional form, in the following we assume a linear 
dependence between the probability of an outcome and the fraction of corresponding monomers 
bound into a complex. These fractions of concentrations are naturally bound between zero and 
one and thus can be linked to probabilities without any need for further renormalization; also, 
such a linear dependence is the simplest form of a monotonic map. Consequently, the fraction of 
target molecules that is bound to the effector, whether such binding is recognized by the guard 
or not, Pt = {[ET] + [RET])/[T]o, quantifies the probability of manipulation of the cellular 
chemistry by the effector. The fraction of the guard protein molecules bound to the effector- 
target complex, Qt = [RET]/[R]q, reflects the efficiency of recognition by the guard protein 
of such manipulations. Likewise, the fraction of the guard protein concentration bound to the 
effector-decoy complex, Qd = [RED]/[R]q, describes the efficiency of recognition by the guard 
protein of the changes induced by the effector in the decoy. The binding between effector and 
decoy by itself does not cause any physiological consequences, hence we do not introduce any 
special notation for the probability of such events. Note that our description takes into account the 
possibility of an "overreaction" by the plant immune system, that is, a possible immune response 
induced by the guard protein when Qt or Qo is large (close to one), while Pt <C 1, i.e. when the 
the actual infection of the cell is rather improbable. 

C. Probabilities of outcomes 

After the pathogen has injected the effector into a cell, three distinct scenarios can develop: 

• Nothing happens to the cell 

• The immune system kills the cell and the infection does not spread 

• The parasite successfully replicates in the cell and subsequently infects other cells 

In terms of the quantities Pt, Qt, and Qd the probability that nothing happens to a cell in the 
presence of the effector is 

Un = {1-Pt){1-Qt){1-Qd). (5) 

This expression means that the effector neither manipulates its target sufficiently to cause the 
infection, nor does it change the target and decoy sufficiently for either change to be recognized 



by the guard protein. Similarly, the probability of infection is given by 

Ui = PT{l-QT)il-QD)- (6) 

Here the effector does manipulate the target sufficiently to cause the infection, yet the changes 
in target and decoy are not recognized by the guard. Finally the probability of activation of the 
immune response and killing the cell, and preventing the spread of infection is 

n/ij = Qt + Qd- QtQd. (7) 

The first two terms on the right hand side describe the probabilities of immune response induced 
either by the effector target or the decoy, and the third term eliminates double counting in the event 
when both the target and decoy are inducing the immune response. These three scenarios form a 
complete set of events, so 

D. Fitness costs 

Now we evaluate fitness costs of these three possible scenarios. First we assume that the cost 
of maintenance of the decoy protein reduces the fitness of an organism by the amount proportional 
to the concentration of this protein, 

n - r [^]o rs^ 

[^ Jo 

Minus signifies the reduction of fitness, and [T]q in the denominator sets the concentration scale. 

An induced immune response implies the death of the attacked plant cell, so the fitness cost of the 

immune response is the cost of replacement of the cell, which we denote by am. 

Finally, the fitness cost of the infection can be expressed through a^ji and a new parameter /i, 

which denotes the number of new cells that become subjected to parasite attack as a consequence 

of successful parasite multiplication in the infected cell. 

0-/ = cTiR + fiijliai + Uirctir). (9) 

The first term on the right-hand side of this self-consistent equation reflects the fitness cost due to 
the death of the infected cell. The second term takes into account that each of fi cells subject to 
parasite attack due to infection of a primary cell either become infected themselves with proba- 
bility 11/ and fitness cost aj, or they successfully thwart the pathogen attack with probability Hir 



and at a cost am due to activation of the immune response. When fiUi > 1, equation ^ has no 

positive solution for a/, which describes the indefinite spread of infection so that all cells in the 

organism die. However, when /ill/ < 1, the organism eventually overcomes the infection with the 

total cost to fitness 

ai = aiR- -— . (10) 

l-Zill/ 

Consequently the total fitness reduction caused by introduction of an effector into a cell is equal 
to the sum of fitness costs, Eqs. (I8I10I) . where the costs of immune response and infection need to 
be multiplied by the probability of the corresponding outcome: 

a = -Cm-— aiRlliR - ct/rII/ IIir + -— . (11) 

[T]o V l-fiUi J 

III. FITNESS OPTIMIZATION 

Once the fitness costs of all three possible outcomes are established, we look at how the organ- 
ism can maximize fitness, or minimize total fitness loss. Since the decoy is assumed to be free 
of other, non-immune roles, its evolution should be the fastest and the least constrained. Thus 
we look for the fitness maximum, varying the parameters related to the decoy which are assumed 
independent: the dissociation constants for binding between the decoy and the effector Ked and 
between R-protein and decoy-effector complex Kred, and the decoy concentration [D]q. Unless 
otherwise specified, all concentrations and dissociation constants are assumed to be constant and 
equal to one, and /x = 5. 

The following general evolutionary pattens were observed when the cost of decoy maintenance 
was varied: 

• The binding between the effector and decoy always tends to be maximized, Ked — > 0, so 
that for both the sink and trigger functions, the decoy absorbs as many effector molecules 
as possible. This seems natural, as it is worthless to maintain any number of decoy copies if 
these copies are not interacting with the effector. Thus in the following we set the effector- 
decoy dissociation constant to Ked = 0.2; the precise value of Ked does not affect our 
conclusion provided that it is noticeably less than all the relevant concentrations, and it 
seems reasonable to assume that evolution of Ked stops when the strong binding limit is 
reached. 
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FIG. 1: Expensive, Cm = 0.3Cc, decoy tends to be less abundant and develop strong binding between 
the effector-decoy complex and the guard protein. In this case the decoy playing the role of an immune 
response trigger yields higher fitness. White area corresponds to death of the organism. 

• For high costs of maintenance of the decoy (roughly for Cm > O.Sa/ij), the fitness optimum 
is reached for smaller decoy concentrations and strong binding (small Kred) between the 
R-protein and effector-decoy complex. Thus, when it is desirable to maintain fewer de- 
coys due to high costs, decoys are predominantly used as effective triggers of the immune 
response. Fig. 1. 

• For low costs of maintenance of the decoy (roughly for Cm < 0.2C(7) the cell produces 
more decoy molecules; decoys simply absorb the effector without triggering the immune 
response, so that Kred is large. Fig. 2. Hence the cheaper and more abundant decoys 
become sinks for effectors, and save more cells from the death inevitably caused by the 
immune reaction. 

• The "survival" curve, which separates the area of the organism death (marked as white 
area at the bottom of all Figures) from the area of survival with a finite fitness cost (grey 
area with contour lines) is practically invariant in all figures and independent on the decoy 
maintenance fitness cost. This means that for small decoy concentration the survival requires 
strong binding between the decoy-effector complex and the guard protein, or, in other words, 
a reliable immune trigger. 
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Kred 

FIG. 2: Less expensive, Cm = O-'iCc, decoy tends to be more abundant and develop weaker binding 
between the effector-decoy complex and the guard protein. In this case the decoy's role of a sink for 
effector provides the higher fitness. 

In a second set of numerical experiments we varied the binding affinity of the guard protein to 
the target-effector complex, Kret- Smaller Kret correspond to a well-tuned immune responses 
that is sensitive to manipulation of the target by the effector, while larger Kret correspond to slow 
and unreliable immune reaction. We observed that: 

• For strong recognizability (Kret = 0.5, Fig. 3), the decoy tends to play the role of a 
trigger, being less abundant (smaller [D]o) and more recognizable for the guard protein 
(smaller Kred)- In this case the probability of inducing immune response is already high 
without a decoy, which makes cell death highly plausible. The additional decoy improves 
the reliability of effector detection by the cell even more, yet the concentration of the decoy 
is not too high due to the finite cost of its maintenance. 

• For weak recognizability (Kret = 2., Fig. 4), the decoy tends to be a more abundant (larger 
[D]o) sink, with the affinity for binding between the decoy-effector complex and the guard 
being smaller (larger Kred)- In this case the cell has more chances to avoid death both from 
the immune response and infection by sequestering the effector through the decoy. The total 
fitness in this case is higher than in the former one. 

Finally, we considered how the severity of the pathogen attack shapes the optimal decoy strat- 
egy. A stronger attack can manifest itself as an elevated concentration of pathogen effectors. In 
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Kred 

FIG. 3: Strong recognizability of effector action on the target by the guard, Kret = 0.5. The decoy tends 
to be less abundant and to make the decoy-effector complex more recognizable to the guard, thus playing a 
role of immune response trigger. 




Kred 

FIG. 4: Weak recognizability of effector action on the target by the guard, Kret = 0.5. The decoy tends 
to be more abundant and to make the decoy-effector complex less recognizable to the guard, thus playing a 
role of a sink for the effector. 

Fig. 5 we show that the increase in concentrations of effector io Eq = 2 leads to stronger recogniz- 
ability (smaller Kred) and smaller concentration of the decoy. A similar effect occurs when the 
concentration of the effector is kept constant, but the effector is more "virulent", so that the affinity 
of binding between the effector and its target increases , Ket = 0.5, Fig. 6. Thus in general, an 
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FIG. 5: The higher concentration of effector, Eq = 2, (left panel) shifts the fitness towards lower decoy 
concentrations Dq and higher recognizability by the guard of the decoy-effector complex (smaller Kred)- 
Right panel shows the control system with £"0 = 1 (same as Fig. 1) , in both panels Cm = 0.2Cc- 

increase in severity of the pathogen attack leads to a switch of the decoy's role from being a sink 
to being trigger (unless the decoy already is a trigger, e.g. due to high maintenance costs). 

IV. DISCUSSION AND CONCLUSIONS 

We developed a model intended to shed some quantitative light on the two basic ways in which 
decoy proteins are thought to function in plant immunity. Decoys, which are most probably du- 
plicated copies of pathogen's effector targets in a plant cell, presumably evolved to play the sole 



role to efficiently mimic effector targets | 



lOh . However, after attracting the pathogen effector that 



invaded a plant cell, the decoys can subsequently function in two ways: In the first role, decoys act 
as sinks for pathogen effector proteins, so that the pathogen attack is thwarted simply by rendering 
the effector proteins useless for changing the hosts cellular environment in favour of the pathogen. 
In the second role, the decoy protein act as triggers for the host immune response, revealing the 
presence of effectors to the guard R proteins that detect the effector-induced changes in the tar- 
gets and decoys, thus effectively making the host immune response more sensitive. So far, these 
two modes of operation were known to occur in different signaling pathways and in different bi- 
ological context: The best studied examples of decoys acting as sinks come from cancer-related 
studies of decoy receptors, which absorb cytokines and other signaling molecules active in the 
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FIG. 6: The stronger binding of effector to target, Ket = 0.5, (left panel) shifts the fitness maximum 
shifting towards lower decoy concentrations Dq and higher recognizability by the guard of the decoy- 
effector complex (smaller Kred)- Right panel shows the control system with Ket = 1 (same as Fig. 1) , 
in both panels Cm = 0.2Cc- 

apoptotic signaling loop [1], and are not directly related to plant immunity. The known examples 
of decoy-triggers come from the rather sophisticated extension of well-established guard-guardee 
plant immune system functioning mechanism dlOD. Our models suggest that the basic functioning 
of these two decoy modes is closely related. In particular, based on our results it can be expected 
that under some conditions, the decoy proteins involved in plant immunity may also play the role 
of simple sinks absorbing pathogen effectors. 

Our quantitative model is based on simple yet realistic biochemical postulates and connects the 
reversible binding between four types of monomer molecules (pathogen's effector, effector target, 
decoy of effector target, and guard protein) to the probabilities of outcomes such as cell death 
and immune reaction, and finally, to the total fitness costs of such outcomes. We assumed that 
the decoy is free to mutate to achieve the highest fitness. Depending on the cost of maintaining 
decoy proteins, on the efficiency of the guard proteins triggering the immune response, and on the 
intensity of the pathogen attack, the optimal decoy strategy may consist of being a sink or of being 
a trigger. When a decoy is costly to maintain, it tends to become a less abundant but a reliable 
immune trigger, while when the cost to express many decoy copies is low, it is more likely to be a 
simple sink. Also, when the immune system functions well without a decoy and the probability of 
detection of pathogen attack is high, the decoy tends to further improve the reliability of triggering 
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the immune reaction rather than to simply absorb effector molecules. Finally, when the severity 
of the pathogen attack increases, either by virtue of higher concentration of effectors, or by more 
efficient binding of effector to its target, the decoy also tends to function as a trigger rather than 
a sink, leaving fewer chances for the infection to pass undetected. Our model also suggests that 
in principle the decoy roles of sink and trigger could continuously adjust to the changing environ- 
ment, both external and intrinsically cellular, through evolution of binding affinities and expression 
levels. 

Presently, we are unaware of the existence of any other quantitative description of the function- 
ing of decoys in the guard-guardee model of plant immunity, whether as a sink or as a trigger. We 
hope that the approach developed here, which links the underlying biochemistry to resulting cost 
of infection and hence to evolutionary considerations, will prove to be useful both for empirical 
and for further theoretical studies of decoy signaling mechanisms in general, and for plant immu- 
nity in particular. For example, it would be nice to have comparative data linking parameters such 
as cost of maintaining decoys and binding affinities in signaling pathways to the functioning of 
decoys as sinks or triggers. 

Our model is based on the assumption that the host species is genetically predisposed for the 
production of decoys, i.e., that the mutation enabling production of decoy proteins has already 
occurred (e.g. through gene duplication). Under these assumptions, the model studies whether 
decoys should be used as sinks or as triggers, and the only scenario in which production of de- 
coy proteins is not favoured at all occurs when such production is too costly. Thus, assuming that 
genetic constraints (e.g. low rate of gene duplication) do not prevent the occurrence of suitable mu- 
tations, decoy mechanisms should be expected to occur in many natural systems. The five decoy 
examples presented in [llOn indicate that this mechanism may have evolved independently in phylo- 
gentically diverse plant species, and hence that decoys may indeed be a wide spread phenomenon. 
The simplicity and potential efficiency of the decoy mechanism suggests that its occurrence may 
not be limited to plants and that it may be used in animal immunity as well, and it remains to be 
seen how ubiquitous and general the use of decoys in immune response systems is. 
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